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Abstract 


Within the context of the restricted problem of tliree bodies, we 
widi to show the effects, caused by varying the mass ratio of the pri- 
maries and the eccentricity of tlieir orbits, upon periodic orbits of 
the infinitesimal mass vdiich are numerical continuations of circular 
orbits in the ordinaiy problem of two bodies. A recursive pouer series 
technique is used to numerically integrate the equations of motion as 
well as the first variational equations in order to generate a t\\^o 
parameter family of periodic orbits and identify the linear stability 
cliaracteristics tliereof . Seven such families are investigated witli 
equally spaced mass ratios from 0.0 to 1.0 and eccentricities of the 
orbits of the primaries in a range 0.0 to 0.6. Stable orbits are 
associated with large distances of the infinitesimal mass from the 
perturbing primar>', nearly circular motion of the primaries and, to 
a slightly lesser extent, small mass ratios of the primaries. On the 
other hand, unstable orbits for the infinitesimal mass are associated 
with small distances fi'om the perturbing primar>', highly elliptic 
orbits of the primaries and large mass ratios. 


I. INTRODUCTION 


In the restricted problem of three bodies one concerns himself 
with tlie motion of a body, Imdng infinitesimal mass, vdiich is sub- 
jected to tlie gravitational influence of two other bodies, called the 
primaries, having finite mass. The assui^ition is made that the body 
of infinitesimal mass exerts no significant gravitational forces upon 
the t»\?o primaries but is only acted upon by them. 

The case where tlie primaries describe circular orbits about their 
center of mass is usually referred to as the Circular Restricted Problem 
of Three Bodies and has been studied by many researchers in great 
detail. Szebehely (1967) provides an excellent text and an extensive 
reference source for most of the work done through 1966 concerning this 
problem, hhere the primaries describe orbits which are ellipses, the 
so-called Elliptic Restricted Problem of Three Bodies, research has 
not nearly been so extensive. See, for example, Szebehely and Giacaglia 
(1964) , Danby (1964a) , Bennett (1965) , Lanzano (1967) , and Broucke (1969) 
Junqueira (Junqueira and Greene, 1970) has been treating tlie restricted 
problem where the motions of the primaries are completely general. 

It is the purpose of tliis paper to present a tivo-parameter survey of 
periodic orbits in the restricted problem of three bodies where the pri- 
maries are constrained to move on circles or ellipses, and the body of 
infinitesimal mass has its motion contained vslthin tlie plane defined by 



those orbits. In particular, we wish to show the effects caused by 
varying the parameters m', the mass ratio of the primaries, and e, the 
eccentricity of the orbit of the primaries, upon periodic orbits of 
the first kind (Poincare, 1892), i.e., periodic orbits of the body of 
infinitesimal mass which are numerical continuations of circular orbits 
in the ordinaiy problem of two bodies. 

Basically, this represents an extension to an earlier paper by 
Shelus and Kumar (1970), which will hereafter be referred to as Paper A. 
Since tlie time of tliat report, mathesnatical computations have been 
improved to eliminate the several inconsistancies which were cited in 
Paper A. Also, the scope of that previous study lias been greatly 
expanded. 


II. H^TIONS OF KJ0T1C»I 


Hie equations of motion for the body of infinitesimal mss in the 
planar elliptic restricted problem of three bodies (Scheibner, 1866; 
Nechvile, 1926; Szebehely, 1967) are usually presented in the following 
fonn: 


where 


x'» - 2y' = (1 + e cos v)’^ (Il-la) 

y" + 2x' = (1 + e cos (Il-lb) 
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^ y^) 


1-u 


[(x4y)‘-+y2]l/2 


[(x-l4uV^+y27l/2 


(II-2) 


Dimensionless units are chosen such that the unit of distance is 
the variable separation between the primaries, the unit of mass is the 
sum of the masses of the primaries (where m is taken to be the mass of 
the less massive of the two primaries, i.e., and the unit of time 

is picked so that the constant of gravitation will be unity. Primes 
denote differentiation vrith respect to the true anomaly, v, of the orbit 
of the primaries, and e is the eccentricity of that orbit. The sub- 
scripts X and y of equations (II-la,b) represent differentiation of B 
witli respect to the indicated variable. Finally, the orthogonal axes 
(x,y) are diosen so that the primaries are located on tlie x-axis when 
V = 0, 


It should be rov-ognized that the equations (n-la,b) are defined 
for a pulsating, nonuniformly rotating coordinate system, where the two 
primaries occupy fixed positions on the x-axis. The origin of coordin- 
ates is at the center of mass of the two primaries and the larger of the 
two masses is to the left of tlie origin. Tliat is to say, ^e coordinates 
of the mass Ini is (-y »0) and those of the mass y is (l*y ,0). Also, 
when e = 0, i.e., when the primaries describe circular orbits, the true 
anomaly of the orbit of the primaries is identical to the dimensionless 
time, and the equations cited reduce to those \diich are craranonly used 
in the planar circular restricted problem of three bodies. Therefore, 
only one set of equations is needed to describe both the circular and 
the elliptic restricted problems. 

At this time one should also note the fact that the equations of 
motion for the elliptic restricted problem do not admit the Jacobian 
integral (Ovenden and Roy, I960) which has been used so extensively in 
studies of the circular problem. 


III. VARIATiam ri^ATIO^S 


In general, the equations of motion (II"la,b) of the previous section 
can alwa>’s be represented by the follcndng system of first order differ- 
ential e«quations 

dx. 

~ f »^2 • * * * i~l,2,...,n. 

The symbol t is used to represent the independent variable instead of 
V, whidi was used 5n Section II. In the elliptic restricted problem, 
the contain only periodic functions of t. 

Let us suppose that the equations (III-l) have periodic solutions: 

X. = j5.(t) = {5. (t+T) (III-2) 

Tlie period T, of course, does not have a uniquely determined value 
since any integral multiple of T is also a period of the solutions. 

By allowing 


x^ = 0^(t) * 


(in-3) 


to be a solution differing only slightly from the known periodic solu- 
tion, we can expand (HI-3) in a Taylor's series, and neglecting the 
squares and higher powers of form the so-called first variational 
equations of the system (Ill-l), i.e.. 


dq 3fj 3f. Sfj 

8fl = 


(III-4) 


Each of the partial derivatives is evaluated along the reference 

j 

solution (I II- 2). 


By the theory of Floquet (Danby, 1964a) , the general solution 
to equations of the form (II 1-4) is 

E, - E s. ,(t) e“j - (III-5) 

1 j-1 

where the S. . are periodic functions of t with period T. The a- 's 

a.t ^ 

are referred to as characteristic exponents and tlie e •' 's are cliar- 
acteristic roots. If all of the characteristic exponents of (III-5) 
are jxire imaginary, the can be expressed as the sums and products 
of purely periodic terms and they will remain finite for all values 
of t. Mien there is a real part in any one of tlie exponents, the 
wHl become unbo’jndcd with t. It should be noted that if all of the 
cliaracteristic exponents were to have negative real parts, the ^ 
would tend to zero as t became large. However, this cannot happen 
for Hamiltonian systems since it can be sham (Pollard, 1966) that 
for those types of systms, if B is a characteristic root, then so 
must -B. 

Following the treatment by Wintner (1947), the variational equa- 
tions (III-4) can be expressed in matrix form, i.e., 

Ci = A(t)^, (I1I-6) 

where A is a inatrLx periodic in t and the and the are each 
column vectors. Dots represent differentiation with respect to the 
independent variable, t. This set of n differential equations will 



be solved when any set of n linearly independent solutions are kncA“TJ. 
If each of these n solutions be a coliBim of a matrix X(t), which is 
called a fimdamental natrix or matrizant (Danby, 1964b) , it is clear 
that 

X(t) = A(t)X(t) (in-7) 

Also, any linear combination of the columns of the fundamental matrix 
is also a solution. That is, Z(t) is also a fundamental matrix if 

Z(t) = X(t)C (III-8) 

where C is seme constant n x n nwtrix sudi that det C 0. 

In particular, let us define one special fundamental matrix, 

i.c. , 

= X(t)[X(tj3'^ (III-9) 

Note thatsXtQ,tQ) is the identity matrix I. Setting X(tg) = I, the 

components of n may be found by ordinaiy^ numerical integration and 

X(t) = n(t^j,t)X(tQ). (III-IO) 

Tlie unique mtrix I2(tQ,tQ+T) is called the mondromy matrix of the 
fundamental matrix X(t) with reference to the given period, T, of A. 

oiT 

The charactcT istic roots, e , of equation (I II -5) are the 
roots of the characteristic equation 


(III-ll) 


- sl| - 0 

i.e. , they are the eigenvalues of the ino: '* X)ny matrix n(tQ,tQ+T). 

” The characteristic exponents are then 

= |log s^ = yllog |s^l + i arg(s^+2kTr)]. (III-12) 

Thte if js^l =1, then log js^| « 0 and yill be pure imaginar>% 

The will therefore be bounded if all of the characteristic roots 
of (I I I- 5) be located on the unit circle in the complex plane. 

In keepir»g with the many other studies which deal with the 
restricted problem of three bodies, we shall use the characteristic 
roots and the diaracteristic exponents, as defined above, to classify 
the periodic orbits which we liave coanputed. Given a periodic orbit 
in the restricted problem, if all four characteristic roots be located 
on the unit circle in the con^lex plane, the orbit will be considered 
to be of the stable t)'pc, i.e., the ^ are bounded, in a linear sense. 

In any other case, the orbit will be considered to be of the unstable 
type. 

Now, returning specifically to the elliptic restricted problem, 
if we let = Ax, ^ = Ay, = Ax', and = Ay' represent the slight 
variations from a known periodic solution, we can obtain, as the first 
variational equations of the system (Il-la,b), the follov\ring expressions: 



(IIM3a) 




= 2(Ay') + (1 + e cos v)'^ 

= -2 (Ax') + (1 + e cos v)’^ 


[jl(Ax) + m(Ay)] 
[m(Ax) + n(Ay)] 


(III-13b) 

(III-13C) 

(III-13d) 
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(x*v)V [(x-h.-l)Vl^^^ (X^-I)V 
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n = 1 ‘ 




^ ^ Vt 
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IV. R]’.sui;j's 


In l!ns .sUiJy :ill i.a'acrical intL-^jrations arc jicrfumc-d usintj '» 
recursive power sorjcs tcclmiquc (Steffciiscn, 19S6; Ilahe, 1903 ; 
l)c]iril and I'ricc, 19(»S; rchlhcrj;. J966; HrourI.c, 3909) as applied to 
botlj llic cciuations of motion, and the first variational 

equations, (III l.vi,b,c,J). It is u-ual to assume that the sum of 
all neglected terns in each scries thus computed Kill be smaller in 
magnitude tlian the first teni! omitted and \<o lia\c obtained vcr>* good 
results by truncating these- series after tt.cnty terms. TliC last few 
tcnius of each series are interrogated to piovidc a criterion for 
integration stejj sir.c. Ibe coimiulaiions were perforr..ed in double- 
precision on a Univac llOS system at the N,\S\ Maiined Spacecraft Center. 

TJic technique which is used to devejope the two paramctCi- survey 
of periodic cibits is a simjilc and straight fonvard one and can be 
outlined in the following rmimcr. .\s a starting point, one selects 
a circular two-body orbit for the l>ody of infinit''siral mass about 
one of tlie primaries, with a i)erioJ equal to some intog’^al ratio of 
the orbital period of tl;c primaries, lissentially, this is a periodic 
orbit in tne circular restricted jiroblem such that -a = 0. 'Jlic initial 
conditions of this first orbit are tlicn used as initial conditions for 
an orbit when n is nc\;, say, 0.01. Tl'.is new orbit will, cf course, 
not be periodic imdor ll*ese new conditions; however, it will not be 
veiy far from being so. 



We use a linear differential correction technique sijnilar to 
that presented by Deprit and Price (1965) to obtain iirqiroved initial 
conditions, A linearly inq>roved solution can, in general, be 
specified by tlic follwing Taylor's expansion: 



(n^-la) 


(IV- lb) 


(IV- Ic) 


(IV-ld) 


We desire i^riaJicity, i.c., we require tliat (x+^)^ » (x+6x)q, 
(y»Ay)j* “ (y+Ay)0, (x+Ax)^ » (x+ax)^, and (y+Ay)^ = (y+^)Q, and 


therefore the follcwving cqpjations can be formed: 



(IV- 2a) 

(IV- 2b) 
(W-2c) 
(IV- 2d) 


where each of ’.he partial derivatives arc obtained through the numer- 
ical integration of the first variational equations. These equations, 
t 4 .y- 2 a,b,c,d), can be solved to provide the corrections which are to 
be applied to tlie previous initial conditions and the process can be 
repeated in an iterative way until the initial conditions of tlie per- 
iodic orbit are reproduced to tlie desired tolerance. For all orbits 
which are presented here, initial conditions after one period are 
reproduced to at least 1 part in tlie tivelfth decimal place. Of course, 
at the same time, we are able to classify this new periodic orbit as 
being "stable" or "unstable" using the characteristic exponents as 
described in the previous section. The entire procedure can tlien be 
repeated where y is now increased to, say, 0.02. In this manner per- 
iodic orbits are obtained for increasing values of the mass ratio of 
the primaries. Remember that all orbits so far described are for the 
circula-r restricted problem, i.e., e = 0. 

At any point in che above procedure, after a periodic orbit has 
been coiqiuted, we can hold the mass ratio constant and, iristead, increase 
the eccentricity of tlie orbit of the primaries from 0.0 to 0.01, say. 

The same iteration procedure is tlien invoked and the result is a periodic 
orbit in tlie elliptic restricted problem with e = 0.01. The eccentricity 
of tlie orbit of the primaries can then be increased to 0.02 and the 
process again repeated. Tlius, starting from a circular orbit in the 
circular restricted problem with y = 0, we can move along increasing 
ecccntricit)' of the orbit of the prijiiarics, or increasing y, to any 
desired value in this step-by-step iterative manner. 


In progressing alor® either e or y , once a nimber of orbits has 
been conpited, it is possible to extrapolate to better approximations 
of the initial conditions of the next orbit by making use of the initial 
conditions of previous orbits. As Broucke (1969) has pointed out, this 
has allowed for a pronounced reduction in the number of iterations which 
are required to arrive at the next orbit. A different two parameter 
’’family" of periodic orbits is produced by starting froa the very begin- 
ning of the procedure witli another circular two -body orbit and then 
preceding in a similar way. 

Hiis report presents the findings for several of these two param- 
eter families of periodic orbits in the restricted problon of three 
bodies.' It is convenient to divide these families into two groups, 
i.e., one group of orbits which enclose both of the primaries and a 
second group of orbits which enclose only one of the primaries (usually 
the more massive of the two). We shall discuss each of these two groups 
separately. 


a- Orbits encircling both primaries 

For tlie cases where tlie orbits of the infinitesimal mass encircle 
both of the primaries, four distinct families have been generated having 
periods of 4 tt, 6tt, 8ti, and lOn, which correspond to periods that are two, 
three, four, and five times, re.pectively, the period of tlie primaries. 
Tlicreforc, these families have been identified by the ratios 2/1, 3/1, 
4/1, and 5/1. Within each family, a particular periodic orbit is 
identified by c, the eccentricity of tlie orbit of the primaries, and y , 


the mass of m 2 , i.e., the primary \%hich is located to the right of 
the origin. The ’’stable-unstable" analysis is sunmarized in Tables 
I, II, III, and IV vtierein the letter S signifies that an orbit is 
"stable", and U, "unstable". Typical orbits are plotted in Figures 
1, 2, 3, and 4; for clarity, these illustrations depict the orbits 
in a baiycentric, inertial coordinate system, not in the rotating- 
pulsating syston vhich has been used for the numerical integrations. 

All of tlie orbits \diich were conqsuted for the t\v^o families iden- 
tified by the ratios 2/1 r.nd 3/1 are of the unstable t)T)e, in the sense 
which has been defined in Section III. These two families evolve from 
circular orbits in the ordinary problem of two bodies which have radii 
of approximately 1.587,,. and 2.080..., respectively. To illustrate 
the nature of these tu’o families of orbits we have deemed it necessar)' 
only to con^jute orbits for mass ratios up to 1/3 0i<0.25) and eccen- 
tricities of the orbits of tlie primaries in a range 0.0 to 0.1. 

The situation is soirowhat more interesting when we consider the 
family of orbits identified by the ratio 4/1. These orbits have 
evolved from a two-body circular orbit having a radius of approximately 
2.519... . Tliey are located farther from the binary system than 

cither of the first t^co families discussed and orbits have been computed 
for the full range of mass ratios, i.e., 0.0 to 1.0. For small values 
of the eccentricity of tlie orbits of tlie primaries, the orbits of the 
infinitesimal mass are linearly stable for all mass ratios. At eccen- 
tricity 0.02, tliosc orbits computed for mass ratios of tlie primaries 


in excess of about 1/3 6i>0.25) prove to be of the unstable type. 

For an eccentricity of the orbit of the primaries of C.03, only the 
orbit for a mass ratio 1/99 (y=0.01) is of the stable type. All 
other orbits which have been computed for this family are linearly 
unstable. 

Finally, we have the family of orbits identified by the ratio 
5/1. The orbits of tliis family have evolved from a circular t\%'o-body 
orbit having a radius of approximately 2.924... . All of these orbits 

are of the stable type. The full range of mass ratios was considered 
and tlie largest value of the eccentricity of the orbit of the primaries 
was 0.09. It is h>'pothesized that these orbits will eventually evolve 
into linearly unstable orbits as the eccentricity is increased further, 
however, some diffiailty has been met in numerically continuing this 
family of orbits aiid further worl: along these lines is required. 

b. Orbits encircling one primary 

Preceding no\i? to the cases where the motion of tlie body of infini- 
tesimal mass encircles only one of the primaries, we are confronted 
with a slightly different situation. Because it is necessary that tlie 
orbital period of the infinitesimal mass be restricted to that which 
is an integral ratio of the orbital period of the primaries, the body 
of infinitesimal mass must circulate about one of the primaries (\diile 
being perturbed by the other) some integral numbers of times during 
one orbital period of the primaries. Tlirce such families have been 
generated having, of course, periods of 2r.. Wc have identified them 
by tlie ratios 1/12, 1/6, and 1/5, i.c., the infinitesimal body will 



conplete twelve, six, and three circuits, respectively, about its pri- 
mary while tlie primaries con 5 >lct one circuit. 

The "stable-unstable" results for these three families are sun- 
marized in Tables V, VI, and VII, where, again, the letter S sig^dfies 
a "stable" orbit and U an "unstable" one. Typical orbits are plotted 
in Figures 5, 6, and 7. These illustrations depict the orbits in an 
inertial syston, havcver, unlike tlie plots for tlie previous four families 
of orbits, the origins here are at the primar>' rather than at the center 
of mass. 

The orbits of the family identified by the ratio 1/12 evolve from 
a circular t\\ 0 -body orbit of approximate radius 0.19078..., such that 
the body of infinitesimal mass encircles its primary exactly twelve 
times during 2ti units of dimensionless time. As is expected, these 
orbits are stable for the entire range of mass ratios when the eccen- 
tricity of the orbit of the primaries is small. Only when this eccen- 
tricity reaches 0.49 do some orbits (for large mass ratios of the pri- 
maries) become unstable. As the eccentricity continues to increase, 
orbits become unstable for smaller mass ratios until, at e = 0.53, all 
mass ratios of the primaries produce unstable orbits for the body of 
infinitesimal mass. 

Continuing to the family of orbits identified by the ratio 1/6, 
we see that these orbits evolve from a two-body circular orbit of 
radius 0.3028S..., with the infinitesimal mass now encircling its 
primar)' exactly six times in 2it units of time. In this case, orbits 
become unstable at smaller orbital eccentricities (e = 0.17) of the 



priwarxes than the previous family discussed, i.e., compare Tables V 
and VI. Also, note that the dependence of instability upon mass ratio 
is far more gradual, such that it is only at e = 0.23 that orbits of 
the infinitesimal mass become unstable for all mass ratios of the 
primaries . 

To conclude this survey we now consider the family of orbits iden- 
tified by the ratio 1/3. The orbits of this family have evolved from a 
circular two-body orbit of radius 0.48074... . Tlie body of infinitesimal 
mass encircles its primar)' only three times in a period of Z-r. units of 
time. All orbits which have been computed prove to be of the unstable 
type, except for a few cases. One group of orbits, for m = 0.23 through 
li - 0.5, are stable for e = 0.0; these orbits become unstable as e is 
increased to 0.02 or greater. Another series of orbits, for e = 0.59 to 
0.65, are stable for m = 0.01 only, and all become unstable when n is 
increased to 0.02 or more. Some difficulties in obtaining convergence 
to orbits for large mass ratios and high eccentricities were encountered 
for this family. These cases are identified in Table VII by tlie symbol "?*’ 



V. Disajssi®^ 


Hie application of the circular restricted problem o£ three bodies 
to Solar Systai dynamics has been n»st iiiqxirtant. Very obvious, of 
course, has been the study of: 1) the equilibrium points and in 

the Sun- Jupiter-Trojan asteroid configuration; 2) the motion of the 
Earth-^toon system as perturbed by the Sun; and 3) more recently, the 
motion of an artificial space vehicle in the Earth -^fc>on neighborhood. 
Somewhat less frequently, see for exanqile liiang and Struve (1956) and 
Ifopal (1959) , the circular restricted problem has been applied to the 
dynamics and the evolution of close binary systens. 

In the study which lias been begun in Paper A and which has been 
continued here, the philosophy for generating periodic orbits in the 
restricted problem of three bodies is perhaps somewhat different than 
that normally assumed in work of this type. We wish to consider this 
survey of orbits in a discussion of the problem concerning the exis- 
tence of low mass Earth-like objects in the vicinity of double stars. 
The components of binary star systems generally describe noncircular 
orbits and their mass ratios exceed those of Jupiter-Sun (^0.001) and 
Moon-Earth ('^.01). It has been pointed out by Iftmiar (1967) that it 
is unlikely that Earth-like objects will survive over long periods of 
time within such binary star systems; a small Earth- like object may 
survive in the vicinity of a doihle star only if it is located close 
to one or to tlie other of the two components, or, if it is located far 
from both stars. 


These remarks must also apply to double star systems in vdiich 
one of the conponents is a ’’dark” ccmpanion. Table VIII gives a 
sunsnary of the observed data for these systems. Note the highly 
elliptic orbits and the large mass ratios of the conponents. The 
dark caipcinions in these systans have been identified witli stars of 
very low mass which have become conpletely degenerate objects, or 
black dwarfs, without going through normal stellar evolution (Kumar, 
1963; 1967). As has been remarked, small Earth-like objects may not 
survive over long periods of time in systems such as Lalande 21185, 
unless they are located in very special positions relative to the two 
primary conponents. 

The results of this surv^ey has at least offered first order sup- 
port for this hypothesis, to the extent that tlie d>Tiamics of small 
masses mthin binary systems can be satisfactorily modeled by the 
restricted problon. Seven families have, so far, been obtained. 

EacJi has evolved from a circular orbit in the ordinary problem of two 
bodies, and eadi orbit encloses one or both primaries. Stable orbits 
tend to be associated with large distances of the infinitesimal mass 
from the perturbing primal^, nearly circular motion of the primaries 
and, to a slightly lesser extent, small mass ratios of the primaries. 
Unstable orbits tend, on tlie other hand, to be associated with small 
distances from the perturbing primary, highly elliptic orbits and 
large mass ratios. 


These results, then, are in keeping with the hypothesis (at 
least in a linear sense) that low mass objects in binajy systems 
will be in ’’stable" orbits when they are located far from both 
primary conponents of the systan or when they are located very 
near to one of the primaries, provided that the eccentricit)’’ of 
the orbit of the primaries is not too large. 
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Linear stability characteristics for orbits of family 
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Linear stability characteristics for orbits of family l/6. 
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TABLE VII 


Linear stability characteristics for orbits of family 1/3 



.0 

.01 

.02 

. 1 

.2 

.3 

.4 

.5 

.6 

.01 

u 

u 

u 

u 

u 

u 

u 

u 

s 

.02 

u 

u 

u 

u 

u 

u 

u 

u 

u 

.03 

u 

u 

II 

u 

u 

u 

u 

u 

u 

.0^1 

u 

y 

u 

u 

u 

u 

u 

u 

u 

.05 

u 

u 

u 

u 

u 

u 

u 

u 

u 

.06 

u 

u 

u 

u 

u 

u 

u 

u 

u 

.07 

u 

u 

u 

u 

u 

u 

u 

u 

u 

.08 

u 

u 

u 

u 

u 

u 

u 

u 

u 

.09 

u 

u 

u 

u 

u 

u 

V 

u 

u 

. 1 

u 

u 

u 

u 

u 

u 

V 

u 

7 

.11 

u 

u 

u 

u 

u 

u 

u 

u 

? 

. 12 

u 

u 

u 

u 

u 

u 

u 

u 

? 

.13 

u 

u 

u 

u 

u 

u 

u 

u 

u 

. i 4 

u 

u 

u 

u 

u 

u 

u 

u 

u 

.15 

u 

u 

u 

u 

u 

u 

u 

u 

u 

.16 

u 

u 

u 

u 

u 

u 

u 

u 

u 

.17 

u 

u 

u 

u 

u 

u 

u 

u 

u 

.18 

u 

u 

u 

u 

u 

u 

u 

u 

u 

. 19 

s 

u 

u 

u 

u 

u 

u 

u 

u 

.2 

s 

u 

u 

u 

u 

u 

u 

u 

u 

.21 

s 

u 

u 

u 

u 

u 

u 

u 

u 

.22 

s 

u 

u 

u 

u 

u 

u 

u 

u 

.23 

s 

u 

u 

u 

u 

u 

u 

u 

u 


s 

u 

u 

u 

u 

u 

u 

u 

y 

.25 

s 

u 

u 

u 

u 

u 

? 

y 

u 

.26 

s 

u 

u 

u 

u 

u 

? 

u 

u 

.27 

s 

u 

u 

u 

u 

u 

? 

u 

u 

.28 

s 

u 

u 

u 

u 

u 

? 

y 

y 

.29 

s 

u 

u 

u 

u 

u 

? 

y 

y 

.3 

s 

u 

u 

u 

u 

? 

? 

y 

y 

: 3 i 

s 

u 

u 

u 

u 

? 

? 

u 

u 

.32 

s 

u 

u 

u 

u 

? 

? 

y 

y 

.33 

s 

u 

u 

u 

u 

? 

? 

y 

y 

. 3^1 

s 

u 

u 

u 

u 

? 

? 

y 

y 

.35 

s 

u 

u 

u 

u 

y 

? 

y 

y 

.36 

s 

u 

u 

u 

u 

y 

? 

y 

u 

.37 

s 

u 

u 

u 

1 ) 

y 

? 

y 

u 

.38 

s 

u 

u 

u 

u 

u 

? 

y 

y 

.39 

s 

y 

u 

u 

u 

y 

? 

y 

u 


s 

u 

u 

u 

u 

y 

? 

u 

u 

1 

s 

V 

u 

u 

u 

u 

? 

u 

y 

J \2 

s 

u 

u 

u 

u 

y 

? 

u 

u 

.'*3 

s 

u 

u 

u 

IT 

u 

? 

? 

y 


s 

u 

u 

u 

u 

u 

? 

? 

y 

. U 5 

s 

u 

u 

u 

u 

y 

u 

? 

y 

.?i6 

s 

u 

u 

u 

u 

u 

y 

? 

y 

.'»7 

s 

u 

u 

u 

u 

u 

u 

? 

u 


s 

u 

IT 

u 

II 

y 

y 

? 

♦ ♦ 

.•’t9 

s 

IJ 

u 

u 

1 ; 

1) 

y 

? 

u 

.5 

s 

u 

u 

u 

u 

II 

y 

7 

y 




table VIII 


Data for systems containing unseen companions 


Star 

1 * 
B* 

e 

Reference 

BD + 20 ^^ 2 k 6 ^ 

1/12 

0 6 

Reuyl (1943) 

Lalande 21185 

1/35 

0.3 

Lippincott (i960) 

Barnard's star 

0.011 

0.75 

van de Kan5» (1969a) 

6 l Cygni 

1/25 

0.5 

Deutsch (i960) 

Ci 23^ 

0.094 

0.9 

Lippincott (I967) 

BD + 6^598 

l/l>l/l 4 

0.6 

Lippincott (I969) 

Tj Cas 

0.044 

- 

van de Kamp (1969b) 

Ci 2547 

0.057 

0.2 

Bieger ( 1964 ) 

* Mass ratio of 

secoidary to 

primary. 


** Eccentricity 

of apparent orbit. 
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Figure S. 
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Figure 7. 


Typical orbits for family 2/1. 
Typical orbits for family 3/1. 
Typical orbits for family 4/1. 
T>'pical orbits for family 5/1. 
Typical orbits for family 1/12. 
T>'pical orbits for family 1/6. 
Typical orbits for family 1/3. 
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